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ABSTRACT 


The  phenomenon  of  cargo  tank  overpressure  during  transfer  operations 
was  analyzed  and  analytical  models  describing  tank  pressure  rise  during 
transfer  operations  validated  with  rcale  model  experiments.  Factors  ex- 
amined affecting  tank  pressure  were  cargo  properties,  transfer  rates,  tank 
characteristics  and  vent  system  design.  Findings  indicate  that  typical 
vent  systems  employed  today  have  adequate  capacity  for  venting  gas  but 
inadequate  capacity  to  vent  liquid  after  the  tank  becomes  liquid  full. 
Furthermore,  it  appears  that  tank  failure  is  inevitable  for  the  case  of 
liquid  overfill  unless  loading  rate  (cu-ft  per  sec)  to  vent  area  (sq--ft) 
ratios  are  kept  below  6 ft/sec.  Currently,  most  transfer  operations 
exceed  this  value  with  tank  failure  expected  less  than  one  minute  after 
the  tank  becomes  liquid  full.  A method  to  evaluate  the  adequacy  of 
existing  cargo  tanks  against  the  overpressure  hazard  is  presented. 
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I.  BACKGROUND  AND  SUMMARY  OF  FINDINGS 

j 

Background  f 

Venting  systems  are  employed  on  marine  vessels  to  relieve  pressure 
differences  which  arise  between  the  cargo  tank  and  the  ambient.  Various 
types  of  venting  systems  are  employed,  depending  on  the  cargo  being  trans- 
ferred. Typical  venting  systems  can  employ  the  use  of  flame  screens,  PV 
valves  and/or  flame  arresters,  as  well  as  other  plumbing  type  fixtures, 
and  may  have  open  or  closed  gauging.  If  the  cargo  transfer  rate  exceeds 
the  vapor  relief  capacity,  or  if  an  accidental  overfilling  of  the  tank 
occurs,  the  resulting  pressure  difference  can,  depending  upon  bulkhead 
structural  design,  result  in  tank  darrage  or  failure. 

An  analytical  study  of  the  overpressure  phenomena  has  been  performed 
by  Arthur  D.  Little,  Inc.  (ADL) ^ . The  analytical  study  evaluates  the 
overpressure  of  any  given  cargo  transfer  operation  based  on  cargo  proper- 
ties. i-ransfer  rate,  tank  structural  characteristics,  and  vent  system 
design.  It  is  the  intent  of  the  investigation  included  herein  to  vali- 
date the  engineering  analysis  of  the  ADL  study  through  the  use  of  scale 
model  experiments  and  to  lend  credence  to  design  and  operational  guide- 
lines to  be  based  on  these  studies. 

Summary  of  Findings 

1.  The  ADL  analytical  formulations  for  the  cases  of  gas  venting  or 
liquid  overfill  during  normal  loading  of  cargo  were  found  to  be  valid 
for  the  range  of  loading  rates,  evaporation  rates,  and  vent  system  re- 
strictions examined. 
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2,  For  normal  vent  systems,  cargo  volatilities,  and  cargo  loading 

rates,  it  appears  that  cargo  can  be  loaded  with  an  adequate  safety 

factor  so  as  not  to  approach  the  average  6 psig  estimated  tank  failure 
1 

pressure . 


3.  An  analytical  formulation  for  offloading  cargo  is  developed  and 
validated  with  scale  model  experiments. 

a.  The  absolute  value  of  the  magnitude  of  negative  pressure 
during  offloading  is  less  than  the  absolute  value  of  the  magnitude  of 
pos  tive  pressure  realized  when  loading  cargo  for  equal  values  of  trans- 
fer rate  and  vent  system  resistance. 

b.  There  is  a possibility  of  the  tank  top  buckling  inward, 
especially  when  the  tank  top  is  cambered,  at  a negative  gage  pressure 
which  is  smaller  in  absolute  magnitude  than  the  6 psig  at  which  the  tank 
would  fail  in  overpressure. 

4.  Simplified  expressions  are  developed  which  can  be  employed  to 
calculate  the  maximum  positive  or  negative  gage  pressure  realized  for 
loading  or  offloading  non-volatile  cargoes. 

5.  For  high  loading  rates  and/or  vent  restrictions,  the  maximum 
tank  top  pressure  changes  rapidly  with  small  changes  in  loading  rate, 
indicating  the  danger  in  employing  these  higher  rates. 

6.  Evaluation  of  the  pressure  buildup  during  liquid  overfill  sug- 
gests that  tank  failure  is  inevitable  unless  very  low  ratios  of  loading 


ate  (ft3/ccc)  to  ver 


area 


} are  employed. 


This  study  found  that 


6 ft/sec  is  the  maximum  permiL  lible  ratio  (based  on  a 4fL/D=10  and  H=8ft) , 
while  ratios  of  30  ft/sec  and  greater  are  commonly  employed. 


7.  The  analytical  model  for  the-  case  of  a blocked  vent  provides  an 


easily  es^iloyed  and  conservative  estimate  Of  the  time  it  will  take  to 
reach  the  tank  failure  pressure. 

8.  Tank  flexibility  and  the  inclusion  of  trapped  air  in  the  tank 
help  to  reduce  the  rate  of  the  pressure  buildup  in  the  tank  during  liquid 
overfill,  but  not  such  that  the  danger  presented  by  liquid  overfill  is 
diminished.  Furthermore,  because  of  the  relatively  short  time  required 
to  reach  tank  failure  (typically  less  than  30  seconds  after  the  tank 
becomes  liquid  full),  it  appears  that  the  overfill  situation  must  be 
prevented  from  occurring. 
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II.  CARGO  TRANSFER  AND  PARAMETERS  AFFECTING  TANK  PRESSURE 

A.  Cargo  Tank  and  Vent  System  Co nfigur atio ns 

Cargo  tanks  on  chemical  tankers  and  barges  are  partitioned  into  indi- 
vidual tanks  for  structural  and  cargo  segregation  purposes.  The  tanks, 
vent  systems  and  pumps  are  matched  to  insure  efficient  loading  and  off- 
loading operations.  Typically,  more  than  one  tank  would  be  loaded  simul- 
taneously, with  groups  of  four  being  most  common. 

Three  types  of  vent  systems  are  the  masthead  system,  the  standpipe 
system  and  the  vapor  recovery  system  (VRS).  Schematic  diagrams  of  all 
three  types  of  recovery  systems  are  presented  in  Figure  1.  In  a typical 
masthead  system,  groups  of  tanks  are  manifolded  to  a contnon  header,  with 
spill  valves,  pressure  vacuum  (PV)  valves  and  flame  arresters  in  the  vent 
system  line.  In  contrast,  a standpipe  system  consists  of  a vertical  pipe 
above  each  tank  to  release  excess  vapor  to  the  atmosphere  and  allows  inde- 
pendent venting  of  each  tank.  The  pipe  is  usually  goosenecked  and 
equipped  with  a PV  relief  valve  and  flame  control  device.  Vapor  recovery 
systems  are  similar  to  the  masthead  system  except  addit  :.onal  piping  is 
added  to  return  the  vented  vapors  to  a shore  disposal  unit.  These  systems 
are  employed  to  eliminate  the  emission  of  pollutant  and  toxic  substances 
to  the  atmosphere. 

All  three  types  of  vent  systems  present  a finite  resistance  to  the 
efflux  of  vapor  when  the  PV  valves  are  open.  Vent  pipes  range  in  diameter 
frc:r,  2 V to  12"  and  are  typically  sized  to  accomodate  the  vapor  displace- 
ment rates.  The  magnitude  of  resistance  to  gas  outflow  is  many  times  ' j 

1 

less  than  that  presented  during  liquid  overfill  when  liquid  cargo  is  j 

flowing  through  the  vent  system.  Further,  it  appears  that  vent  systems  j 

] 

1 
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designed  for  gas  venting  cannot  handle  the  liquid  overfill  situation  un- 
less pumping  rates  are  extremely  low. 

A convenient  index  of  the  resistance  to  either  gas  or  liquid  flow  is 
the  effective  length  to  diameter  ratio  (L/D)  of  the  vent  system.  The  1,/D 
of  a vent  system  can  1 at  by  summing  the  D/D  ratios  cf  all  the 

pipe  lengths  employed  ii  -ystem  and  then  adding  in  the  appropriate 

handbook  L/D  values  for  the  various  pipe  fittings,  valves,  flame  arrestors 
or  other  items  which  comprise  the  overall  L/D  of  the  vent  system.  When 
the  L/D  ratio  is  multiplied  by  the  friction  factor  4f  (accounting  for  pipe 
wall  roughness  and  flow  rate)  the  overall  frictional  resistance  of  the 
vent  system  is  established  and  can  be  employed  in  pressure  drop  calcula- 
tions . 

Due  to  the  long  lengths  of  pipe  and  the  number  of  valves  and  arrestors 
ertployed,  vapor  recovery  systems  have  the  highest  L/D  ratios.  Length 
to  diameter  ratios  on  the  order  of  1000  would  not  be  unrealistic  for  a 
typical  VRS  configuration,  while  a typical  manifolded  or  standpipe  vent 
system  would  have  an  I./D  ratio  on  the  order  of  several  hundred. 

B,  Factors  Contributing  To  Pressure  Rise 

In  general,  typical  vent  systems  can  adequately  discharge  the  vapor/ 
air  mixture  developed  inside  the  tank  without  creating  excessive  tank 
pressures  (usually  less  than  a few  psig) . However,  in  the  case  of  an 
accidental  overfill,  or  a stuck  valve,  vent  pipe  blockage,  or  excessive 
loading  or  evaporation  rates,  excessive  tank  pressures  can  occur.  The 
: various  factors  important  to  pressure  rise  during  normal  cargo  transfer 

are  depicted  in  Figure  2.  As  can  be  seen  from  Figure  2,  the  cargo  being 
loaded  is  evaporating  as  it  is  being  loaded  and  the  entering  liquid 
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FIGURE  2 SCHEMATIC  DIAGRAM  SHOWING  FACTORS  IMPORTANT  TO  PRESSURE  RISE 

(FROM  REFERENCE  1) 


displaces  the  vapor/air  mixture  through  the  vent  system.  The  finite  re- 
sistance offered  by  the  vent  system  to  the  efflux  of  the  gas  mixture  cre- 
ates the  resulting  pressure  rise  inside  the  tank.  During  liquid  overfill 
it  is  the  resistance  of  the  vent  system  to  liquid  outflow  that  results 
in  the  pressure  rise  inside  the  tank. 

As  was  stated  earlier,  overpressure  can  occur  by  inadvertent  overfill, 
excessive  loading  rates,  or  inadequate  vent  system  capacity.  Some  of  the 
parameter  values  which  constitute  normal  and  extreme  situations  for  cargo 
transfer  are  given  in  Table  1 (From  Wilson  and  Raj , reference  (1)). 

Further  use  of  the  terms  normal  or  typical  will  reflect  the  normal  values. 

Wilson  and  Raj1  also  made  estimates  of  the  minimum  internal  tank 
pressure  loadings  required  to  initiate  failure  of  cargo  tanks  for  three 
representative  vessel  designs:  an  offshore  barge,  an  inland  barge,  and 
a large  tankship.  Cn  the  basis  of  the  analyses  of  the  three  vessels,  it 
has  been  calculated  that  an  average  internal  pressure  level  in  the  cargo 
tanks  of  about  6 psig  will  be  sufficient  to  initiate  failure  of  the  tank 
structure.  This  nominal  pressure  level  represents  the  average  of  about  8 
psig  for  the  tankship,  6 psig  for  the  offshore  barge,  and  about  4 psig 
for  the  inland  barge. 


TABLE  1 

PARAMETER  VALUES  FOR  NORMAL  AND  EXTREME  TRANSFER  OF  CARGO 
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III.  EXPERIMENTAL  MODEL  FOR  CARGO  TANK/VENT  SYSTEM 


. Analytical  Models  of  Tank  Pressure  pise 


Hie  rationale  behind  the  development  of  the  experimental  models 
can  best  be  seen  by  first  considering  the  analytical  model  of  the  physical 
system  and  the  factors  important  to  pressure  rise.  Details  of  the  deriva- 
tion of  the  analytical  model  can  be  found  in  Wilson  and  Raj^.  Briefly, 
a mass  balance  equation  was  written  for  the  vapor  mass  in  the  tank  as  a 


function  of  the  outflow  rate  and  the  cargo  evaporation  rate: 


Rate  of  change 
of  air/vapor  mass 
in  the  tank 


Mass  rate 
of 

venting 


Mass  added 
from 

evaporation 


Equation  (1)  is  the  basic  equation  describing  the  pressure  rise  in- 
side the  tank.  Then,  by  deriving  an  expression  for  the  evaporation  rate 
and  taking  into  account  the  purping  rate,  vent  system  hydraulics  and  the 
fluid  characteristics,  the  following  equation  was  arrived  at  describing 


the  rate  of  pressure  rise  in  the  tank: 


^Si>.  p/p  .JiSiiiV1  - ,[(p/^’2  ~ l]  ^ . ( 

fill  fill  a ^ t j (Q/A)  [21n(p/pa)+4f  L/d] 


where : 


t = tank  fill  time  (tank  volume/ loading  rate) 
t = time 


p = tank  top  pressure 


= atmospheric  pressure 

K = ratio  of  evaporation  rate  to  loading  rate,  taken  at  t^^ 
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R = gas  constant 

T = temperature 

Q = loading  rate 

A = vent  pipe  cross  sectional  area 

f = fanning  friction  factor 

L/D  = length  to  diameter  ratio  of  the  vent  system 
Equation  (2)  is  subject  to  the  initial  condition  p/p  = 1 at  t = 0 and 
can  be  numerically  integrated.  The  key  independent  parameters  in  equa- 
tion (2)  are  Q/A,  K,  4f  L/D  and  t/t  . 

Formulation  of  the  analytical  model  for  the  case  of  liquid  overfill 
was  approached  .in  the  same  manner  as  that  for  the  normal  transfer  case. 
The  physical  factors  for  modeling  the  tank  pressure  rise  during  liquid 
overfill  are  shown  in  Figure  3.  When  the  tank  becomes  liquid  full,  con- 
tinued pumping  begins  to  displace  liquid  into  the  vent  system.  The 
frictional  pressure  drop  from  the  vent  system  in  addition  to  the  liquid 
head  it  imposes,  creates  the  pressure  buildup  in  the  cargo  tank.  The 
tank/vent  system  attenpts  to  relieve  the  pressure  drop  by  compression 
of  the  liquid  and  expansion  of  the  tank  walls. 

To  analytically  model  liquid  overfill,  a mass  balance  for  the  tank- 
liquid  system  is  first  written: 


dM 

dt 

- 

M. 

1 

f>m 

(3) 

Rate  of  accumulation 

/ ■s 

Inlet  mass 

"outflow  ratc^ 

of  mass 

flow  rate 

of  liquid 

within  the  tank 

of  liquid 

from  the  tank 

where  p is  fluid  density,  A is  the  vent  cross  section  area,  and  U 
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Cargo  Fills  Tank 


t = 1 - 10  sec.  t = 20  - 200  sec. 

Cargo  Enters  Vent  Cargo  Fills 

Vant  System 


FIGURE  3 SCHEMATIC  OF  THE  PHYSICAL  SYSTEM  FOR 
MODELING  THE  TANK  PRESSURE  RISE 
DURING  OVERFILL  (Wilson  and  Raj  (1)  ) 


is  the  mean  velocity  of  liquid  in  the  vent  pipe.  Now,  the  hydraulics 
of  the  vent  system,  the  compressibility  of  the  liquid,  the  cargo  loading 
rate,  and  the  tank  geometric  and  structural  characteristics  are  employed 
to  arrive  at  the  following  governing  equation  for  liquid  overfill j 


dp  _ 1 

dt  ~ (B  +x") 


where : 


(Q/A)  - u 


. -c 

(V^/A)  + (Q/A)  t - J U dt 


(4) 


1 dV 

B = coefficient  of  volume  expansion  of  the  tank  by  pressure  {=—  — ) 

VT  dP 

= 1 d£ 


compressibility  of  liquid 


V = tank  volume 
T 


/o  dp 


and  all  other  terms  are  as  previously  defined.  This  basic  equation  for 
overfill  employed  the  following  assumptions  for  analytical  solutions. 

(1)  Constant  mass  inflow  rate  M 


l 

(2)  Constant  values  of  B andX 

(3)  Fully  turbulent  pipe  flow 

Equation  (4)  is  rather  cumbersome  and  an  iterative  procedure  must  be 
employed  to  obtain  p and  U as  functions  of  time.  The  case  of  a blocked 
vent  provides  a conservative  estimate  of  the  time  required  to  reach  maxi- 
mum allowable  pressure  for  normal  and  high  loading  rates  during  overfill. 
Assuming  U = 0 (vent  blocked)  equation  (4)  becomes: 


P 


1 

(B  +J(> 


In 


(5) 


Saic  rctt-Os? 


at  which  overfill  can  be  sustained  can  be  found  by  assuuiiny 


a steady  flow  condition  (Q/A=U)for  the  vent  system.  Then  the  tank  pressure 
can  be  related  to  the  mass  outflow  rate  and  solved  for  the  liquid  efflux 


velocity  U: 


Q/A  » U - -q 


L/D+gJ 


where  H is  the  maximum  height  of  the  vent  system  above  the  tank  top.  By 
substituting  the  maximum  allowable  gage  pressure  (p-p  ) that  a tank  can 

Si 

withstand  into  equation  (6) , one  would  be  able  to  arrive  at  the  maximum 
allowable  ratio  of  loading  rate  to  vent  cross  sectional  area  (Q/A)  that 
a tank  is  capable  of  receiving  for  an  overfill  situation  without  failing 


(2£kE 


rimental  Model  for  Normal  Transfer 


Physically,  the  test  fixture  consisted  of  puirps  and  flowmeters  for 
air  and  water,  a model  tank  and  vent  system,  and  various  pressure,  strain 
and  deflection  transducers  arranged  in  the  manner  shown  in  Figure  4.  The 
model  tank  is  a 55  gallon  drum  and  the  model  vent  system  is  3/16  inch 
inner  diameter  high  pressure  tubing  mounted  in  one  of  the  drum  bungs.  It 
should  be  noted  that  the  analytical  model  for  normal  transfer  is  in  no 
way  dependent  on  tank  geometric  or  structural  characteristics.  The  choice 
of  a 55  gallon  drum  for  the  model  tank  was  based  on  economy,  convenience, 
and  on  the  drum  coefficient  of  volume  expansion  B,  which  governs  the  rate 
of  pressure  rise  for  liquid  overfill.  The  value  of  B for  the  model  tank 
is  .000388  psi  1 which  is  within  the  range  of  B values  for  typical  +*r>v-g 
,(.001  to. 00001  psi  ) . 

Examination  of  equation  (2)  shows  the  pressure  rise  inside  the  tank 
to  be  a function  of  the  quantities  Q/A,  4fL/D,  t/tf^  RT  and  K (for  the 
case  of  K/0) . Thus  for  any  experimental  model  study,  if  these  quantities 
are  held  the  same  in  the  model  as  in  the  prototype,  the  resulting  pressures 


The  modeling  of  each  of  the  above  parameters  was  accomplished  as  follows: 

(1)  Gas  Properties (RT)  - The  analytical  model  assumes  that  the  vapor- 
air  mixture  being  vented  has  the  same  molecular  weight,  specific  heat  and 
viscosity  as  that  of  pure  air.  Thus,  an  experimental  model  venting  air 

at  the  same  rate  at  which  the  vapor-air  mixture  is  vented  in  the  analytical 
model  would  suffice  for  the  validation  of  the  analytical  model. 

(2)  Vent  Pipe  Efflux  Velocity  (Q/A)  - This  term  is  expressed  as  the 
ratio  of  the  tank  loading  rate  to  vent  cross  sectional  area,  and  serves 
as  a convenient  way  to  express  a nominal  efflux  velocity  based  on  opera- 
tional and  design  guidelines.  Thus,  by  manipulating  this  ratio  (Q/A),  the 
effects  of  different  loading  rates  on  tank  pressure  can  be  examined,  with 
the  resulting  experimental  model  tank  pressures  being  the  same  as  those 
that  would  be  expected  for  a real  life  situation  with  the  same  0/A  ratio. 

(3)  Vent  System  Frictional  Flow  Resistance  - This  parameter  is  taken 
into  account  in  the  term  4f  I/D.  Thus,  for  a given  vent  cross  section 
area  (defines  D)  and  knowing  the  relative  roughness  of  the  model  vent  pipe 
(establishes  4f) , one  can  then  employ  different  lengths  of  pipe  in  the 
tank/vent  system  model  to  evaluate  the  quantity  4f  L/D.  Establishing 
the  model  vent  pipe  relative  roughness  and  relating  it  to  the  fanning 
friction  factor  ( 4 f ) was  accomplished  in  a seperate  determination  and  is 
detailed  in  Appendix  A.  Obviously,  the  flow  in  the  vent  system  must  be 

in  the  same  regime  (e.g. , turbulent)  for  both  the  experimental  and 
analytical  models.  Also,  subsonic  flow  must  be  maintained. 


(4)  Cargo  Evaporation  Rate  - Cargo  evaporation  during  the  loading 
of  volatile  cargo  was  simulated  by  pumping  pir  into  the  model  tank  at  a 
mass  rate  described  by  the  analytical  model  for  cargo  evaporation 
(equation  (5),  Reference  1).  The  mathematical  expression  governing 
evaporation  in  the  analytical  model  is  given  as: 

"vap  " K/>a  « (7) 

and  is  illustrated  in  r-'.gure  5 for  a loading  rate  Q/A  = 100  ft/sec  and 
evaporation  rate  to  loading  rate  ratio  of  K = .2.  The  term  K,  the  ratio 
of  the  cargo  • Dlumetric  evaporation  rate  (M vap^VaP  t0  loadin9  rate  Q 1S 
simply  a means  of  relating  the  evaporation  rate  to  the  loading  rate  at  a 
particular  time  (t  = tfi^)  since  the  evaporation  rate  changes  with  time. 

The  following  assumptions  apply  to  the  analytical  evaporation  model 
and  would  likewise  apply  to  any  experimental  results  where  evaporation  is 
considered: 

(1)  The  density  of  pure  vapor  is  the  same  as  that  for  pure  air 
(/°a  "/W* 

(2)  Cargo  loading  temperature  is  close  to  ambient 

(3)  The  liquid  surface  stays  at  the  cargo  loading  temperature  and 
has  a constant  area  equal  to  the  area  of  the  tank  floor. 

Also  shown  in  Figure  5 is  the  experimental  simulation  of  the  analytical 
evaporation  model  for  the  same  cases.  The  problem  of  M,  going  to  in- 
finity when  t=0  was  approximated  by  setting  the  air  flow  rate  to  the 
maximum  amount  the  air  flowmeter  could  accurately  monitor  until  an 
value  was  reached  which  corresponded  to  an  amount  predicted  by  the  analy- 
tical model.  The  analytical  curve  was  then  approximated  in  the  stepwise 
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manner  shown  m Figure  5.  Additionally , the  pumping  of  air  to  simulate 
evaporation  is  consistent  with  the  assumption  in  the  analytical  model 
that  pure  air  is  being  vented  through  the  vent  system. 

Then,  with  the  above  quantities  modeled  as  described,  one  would  ob 
tain  tank  pressures  in  the  model  which  would  be  the  same  as  prototype 
pressures  that  would  be  realized  in  a prototype  with  the  same  values  of 
the  independent  parameters  4fl/D,  Q/A, and  K,  at  equal  values  of  t/t^_  ^ 
C.  Experimental  Models  for  Liquid  Overfill  and  Blocked  Vent 


The  model  for  liquid  overfill  is  physically  the  same  as  that  for  the 
case  of  normal  transfer.  One  will  note  from  equation  (4)  that  the  rate 
of  pressure  rise  for  liquid  overfill  is  a function  of  loading  rate  Q, 
tank  volume  V , liquid  corqaressibility  vent  area  A,  the  tank's  co- 
efficient of  volume  expansion  B and  the  velocity  of  liquid  in  the  vent 
system  u,  which  is  itself  a function  of  p and  the  physical  configuration 
of  the  vent  system  (i.e.  4f  L/D  and  H,  the  vent  syot-m  resistance  and 
height  respectively) . It  can  be  shown  that  for  mos  c.  ks,  the  bulk 
modulus  B is  much  larger  (less  stiff)  than  the  compret  ability  of  the 
liquid  (B»X)  , therefore  B is  the  more  influential  of  these  two  terms  in 
governing  the  rate  of  pressure  rise.  Then,  if  one  were  to  select  a model 
tank  with  a value  of  B similar  to  that  for  a typical  cargo  tank,  in  addi- 
tion to  ratios  of  loading  rate  to  tank  volume  (Q/V^) , loading  rate  to 
vent  area  (Q/A) , 4fL/D  and  H which  correspond  to  typical  cargo  operations 
one  would  obtain  a model  pressure-time  history  which  would  correspond 
exactly  to  a typical  cargo  tank  pressure-time  history  with  the  same 
values  of  Q/V  , Q/A,  4f L/D,  H and  B.  The  terms  4fL/D  and  H,  while  not 
in  equation  (4) , along  with  p determine  the  rate  of  liquid  flow  in  the 


vent  o.  Thus,  any  comparison  between  experimental  and  analytical  results 
for  liquid  overfill,  should  take  into  account  the  vent  system  parameters 
4f h/D  and  H to  have  a meaningful  comparison.  Additionally,  it  is  worth 
noting  that  the  pressure-time  relationship  in  model  and  real  life  for 
liquid  overfill  is  based  on  the  time  after  the  tank  becomes  liquid  full, 
and  not  the  non-dimensionalized  fraction  of  total  fill  time  as  in  the 
case  of  normal  transfer. 

Thus,  a model  tank  was  needed  which  possessed  structural  characteris- 
tics resulting  in  a modulus  of  pressure/ volume  expansion  within  the  range 
of  typical  moduli  one  would  find  for  an  actual  cargo  tank.  After  in- 
vestigating various  possibilities,  it  was  determined  ttot  a 55  gallon  drum 
possessed  these  characteristics.  Since  95%  of  the  volume  change  of  the 
drum  from  pressure  occurs  in  the  drum's  top  and  bottom,  it  was  imperative 
that  the  drum  top  and  bottom  plate  deflections  be  measured  accurately. 

To  acconplish  this,  the  Moire  Contour  Sum-Contour  Difference  Method 
(References  2 ar.d  3)  was  used  to  measure  the  plate  deflections,  and  was 
backed  up  by  linear  potentiometers  mounted  on  the  drum  top  and  bottom. 

The  Moire  technique  determines  changes  in  surface  elevation  and  hence 
volume,  by  observing  contour-like  patterns  created  by  the  mode  inter- 
ference of  a grid  with  its  shadow  cast  onto  a surface  and  comparing  the 
interference  patterns  for  different  deformed  states.  The  Moire  method 
was  employed  both  statically  (to  calibrate)  and  dynamically  (during  over- 
fill tests)  , for  determining  the  voluuie  changes  in  the  tank  with  pressure. 

✓ 

Details  on  the  technique  and  results  from  the  Moire  evaluations  can  be 
found  in  Appendix  B. 
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Simulation  of  a blocked  vent  was  accomplished  by  installing  a ball 
valve  in  the  vent  pipe  opening  and  closing  it  after  overfill  had  begun. 
As  can  be  observed  from  the  analytical  formulation  for  a blocked  vent, 
equation  (5)  , the  rate  of  pressure  rise  is  a function  of  the  tank  bulk 
modulus  (B),  liquid  compressibility  (d<)  and  the  ratio  of  loading  rate  to 
tank  volume  (Q/V, ).  By  modeling  these  parameters  in  the  same  manner  as 
was  done  for  liquid  overfill,  pressure-time  histories  the  sane  as  those 
ejected  from  a prototype  with  the  same  parameters  would  be  obtained. 
Additionally,  since  the  model  tank  was  rated  for  a maximum  internal 
pressure  of  40  psig,  all  liquid  overfill  evaluations  were  terminated  at 
an  internal  model  tank  pressure  of  25  psig. 
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IV.  PRESSURE  HISTORIES  LOADING  CARGO 


A.  Loading  Non-Volatile  Cargo 

The  effect  of  loading  a non-volatile  cargo  was  simulated  by  pimping 
water  into  the  test  tank  at  a rate  to  produce  the  desired  Q/A  ratio  and 
monitoring  the  tank  pressure  buildup.  The  effects  of  vent  restriction 
(4f  L/D  = 10,  20,  35)  and  loading  rate  <Q/A  = 30,  50,  100,  150  ft/sec) 
with  no  cargo  evaporation  (K  = 0)  were  examined.  Figures  6 and  7 show 
the  pressure  rise  under  the  conditions  of  Q/A  = 100  ft/sec  and  4f  L/D 
= 10  and  4f  L/D  = 20  respectively.  The  tank  top  pressure  (psig)  is 
plotted  against  the  fraction  of  total  fill  time  in  these  figures.  Also 
included  in  Figure  6 is  the  pressure-time  history  as  predicted  by  the 
analytical  model.  One  can  observe  from  both  figures  the  good  repeata- 
bility the  experimental  results  exhibited  from  run  to  run  and  also  note 
the  good  correlation  with  the  analytical  results.  There  is  a slight  dis- 
crepancy between  experimental  and  analytical  values  in  the  early  stages  of 
loading,  but  the  important  point  to  note  in  Figure  6 is  that  both  analyti- 

m 

cal  and  experimental  models  predict  the  same  maximum  pressure  that  will 
be  realized  for  the  loading  operation. 

Equation  (2),  which  describes  the  rate  of  pressure  rise  as  a function 
of  time  must  be  numerically  integrated  and  is  therefore  not  well  suited 
for  practical  use.  A similar,  but  less  complicated  solution  to  equation 
(2)  when  considering  non-volatile  cargoes  can  be  arrived  at  by  letting 
t = and  K = 0 to  arrive  at  an  expression  which  describes  the  tank 

top  pressure  for  loading  cargo  with  no  cargo  evaporation  as  was  done  to 
arrive  at  equation  (8) . 
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K TOP 


•r'*i 


P RT 

Pa  RT  - <Q/A)‘i[21nCp/pa)+4fL/D] 


Equation  (8)  predicts  the  tank  top  pressure  at  t = which 

experimental  and  analytical  pressure-time  histories  indicate  is  the 

maximum  pressure  that  will  be  attained  for  a loading  operation  with  no 

cargo  evaporation.  Equation  (8)  cannot  be  solved  directly,  but  must  be 

iterated  to  arrive  at  a solution.  The  iterative  process  is  hastened  by 

first  solving  equation  (8)  for  p/pa  without  the  term  2 ln(p/p  ) in  the 

denominator  and  then  using  the  value  of  p/p  in  the  complete  expression 

a 

to  generate  a new  p/p  . Convergence  is  rapid  since  2 In (p/p  ) is  small 

a a 

compared  to  4f  L/D.  Table  2 gives  the  experimen tally  arrived  at  maximum 
pressures  for  the  loading  rates  and  vent  restrictions  examined.  Also 
included  in  Table  2 are  the  maximum  predicted  pressures  using  equation 
(8) . The  correlation  between  analytical  and  experimental  results  is 
considered  good  and  certainly  within  the  bounds  of  experimental  error. 

Figure  8 illustrates  the  relationship  between  the  maximum  absolute 
gage  pressure  that  would  he  realized  for  a particular  transfer  rate  (Q/A) 
and  vent  system  (4f  L/D)  as  predicted  by  equation  (8).  From  Figure  8 it 
is  evident  that  with  higher  loading  rates  and  vent  system  frictional 
resistances,  one  need  only  vary  the  loading  rate  slightly  to  produce  large 
changes  in  the  maximum  pressure  attained,  implying  pumping  rates  must  be 


accurately  monitored  or  that  these  rates  must  be  avoided  to  diminish  the 
risk  of  accidental  overpressurization.  Further,  4f  L/D  is  subject  to 
change  (increase)  as  corrosion  develops  in  the  vent  system.  Therefore, 


TABLE  2 


ANALYTICAL  AND  EXPERIMENTAL  MAXIMUM  PRESSURES 
LOADING  NON-VOLATILE  CARGO 


P 

rax 


(psig) 


Q 

A 

( ft/sec) 

NOMINAL  4f^ 

* (ACTUAL) 

10 

* (ACTUAL) 

20 

* (ACTUAL) 

35 

Experimental 

Analytical 

ao_  (12-6) 

.09 

.14  (25*2) 

.19 

,17.  (44-2) 

. 34 

Experimental 

Analytical 

.25  (11*4) 

~ 24 

.52  (22*8) 

T 49 

.75  (41’0) 

.92 

„ Experimental 

Analytical 

.90  (10-0) 

790  " 

1.90  (I9‘5) 

1797" 

3.80  (34’0) 

4 .“o  2" 

Experimental 

Analytical 

1.95  (9*2) 

2.12 



J.57.  (18'1) 
5.44 

** 

* Actual  4f—  value  used  experimentally  and  for  calculations 
(see  Appendix  A) . 

**  Loading  rate  and  vent  restriction  combination  produce 
choked  flow. 
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safe  loading  rates  for  a given  system  may  become  unsafe  in  time  as 
corrosion  may  increase  the  value  of  4£  t/D. 


From  the  analyses  and  results,  the  following  conclusions  were  reached 
for  loading  non-volatile  cargoes: 

(1)  Maximum  tank  top  pressures  resulting  from  loading  cargo  without 
evaporation  are  sufficiently  low  so  as  to  provide  an  adequate  safety 
factor  for  currently  employed  loading  rates  and  tank  structures  $ = 80ft/ 

A 

sec,  4f  — = 20) . 

(2)  Maximum  tank  top  pressures  occur  at  t = t^  for  loading  cargo 
with  no  cargo  evaporation,  and  can  be  calculated  using  the  simplified  ex- 
pression, equation  (8). 

B.  Loading  Volatile  Cargo 

loading  cargo  with  the  cargo  evaporating  into  the  tank  was  simulated 
by  simultaneously  pumping  air  and  water  into  the  model  tank  as  described 
previously.  The  c l lowing  parameter  values  were  examined  to  assess  the 
effects  of  cargo  evaporation: 

Vent  restriction  = 4f  L/D  = 10,  20 

Loadi.ng  Rate  = . _ 

Vent  Area  ^ 

Evaporation  Rate  = R ^ 5 

Loading  Rate 

In  general,  the  effect  of  cargo  evaporation  is  to  increase  the  maximum 
tank  pressure  realized  for  a given  vent  restriction  and  loading  rate. 
Comparing  the  experimental  results  presented  in  Figures  9 and  10  shows  the 
effects  of  increasing  evaporation  and  loading  rates  on  the  pressure-time 
history.  Equation  (2)  shows  that  dp/dt  is  greater  for  evaporation  than  non' 
evaporation  cases.  Figure  11  compares  pressure-time  histories  generated 


EVAPORATION  RATE 
K “ LOADING  RATE 

4f^  = 10 

D 


0.2 


50 


ft/sec 


30  ft/sec 


(OlSd) 


by  the  analytical  and  experimental  models  when  considering  cargo  evapora- 
tion. The  results  shown  in  Figure  11  are  for  typical  evaporation  and 
loading  rates  and  are  indicative  of  the  repeatability  of  the  experimental 
results  and  the  correlation  attained  with  analytical  data. 

The  results  of  the  experimental  evaluations  of  loading  cargo  with 
cargo  evaporation  are  summarized  in  Table  3.  Table  3 gives  the  maximum 
tank  top  pressures  attained  loading  cargo  at  various  transfer  rates, 
evaporation  rates  and  vent  system  restrictions.  One  will  note  that  the 
maximum  tank  top  pressures  attained  are  below  the  tank  failure  pressures 
calculated  for  typical  cargo  tanks  with  the  exception  of  the  case  where 
Q/A  = 100  ft/sec,  K * .5  and  4f  L/D  = 20,  which  combines  high  evaporation 
and  loading  rates.  Thus,  one  can  conclude  that  the  danger  of  tank  over- 
pressure during  normal  cargo  transfer  with  typical  vent  restrictions  and 
transfer  and  evaporation  rates  is  small  and  will  only  occur  if  one  or 
more  of  these  parameters  becomes  excessive. 

C.  Summary  of  Results  loading  Cargo 

The  experimental  results  obtained  for  loading  cargo  are  summarized 
in  Figures  12  and  13  as  plots  of  maximum  tank  top  pressure  vs.  cargo 
trasfer  rate  for  specific  vent  restrictions  and  evaporation  rates.  Also 
included  in  Figure  12  are  the  analytical  results  for  the  parameter  values 
corresponding  to  the  experimental  results,  again  demonstrating  the  good 
correlation  between  the  two.  The  regime  of  expected  tank  failure  is  also 
indicated  in  each  figure  and  one  will  observe  that  - cakes  a combination 
of  high  loading  and  evaporation  rates  to  approach  a tank  pressure  that  is 
considered  unsafe. 
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•,  O EXPERIMENTAL 
ANALYTICAL (EQN  12,  REF  1) 


OltH 


V.  EVALUATION  OF  TANK  PRESSURE  UNLOADING 


CARGO 


A.  Analytical  Formulation 

The  following  formulation  describing  tank  top  pressure  while  un- 
loading cargo  was  arrived  at  in  a manner  similar  to  that  employed  by 
Wilson  and  Raj1  in  developing  an  analytical  model  describing  tank 
pressure  as  a function  of  time  when  loading  cargo.  The  effect  of  cargo 
evaporation  on  the  rate  of  pressure  rise  was  not  considered  since  any 
cargo  evaporation  would  tend  to  reduce  the  magnitude  of  negative  gage 
pressure  created  inside  the  tank. 

A mass  balance  on  the  vapor/air  mixture  in  the  tank  is  first 
written: 


dM 

dt 


M 


'Rate  of  change'  Mass  rate 

of  air  mass  of 

in  the  tank  venting 


(9) 


The  left  hand  side  is  developed  as  follows.  The  mass  in  the  tank  for 


a compressible  gas  can  be  expressed  as : 


" V*  ‘ (£Vt)  <10> 

where : 

/>  - density  of  vapor  air  mixture 

Pa  = density  of  gas  or  air  at  initial  condition 

v = volume  of  vapor  space  at  any  instant  of  time 

V = total  tank  volume 
T 

Letting  be  the  initial  tank  volume  not  occupied  by  liquid  and  tg 

be  the  characteristic  unloading  time  (t  = (V  - V.)/0),  then: 

e T 1 
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(11) 


v_ 


V.  + Qt 


Additionally,  for  an  isothermal  gas,  AA  = P/P&  and  equation  (10) 
can  be  expressed  as: 


M =/>  (V  + Qt  ) (P/P.) 

• a A e a 


Vi  + Qt 


V.  + Qt 
i e 


(12) 


The  venting  rate  is  now  related  to  the  pressure  rise.  Engploying 

4 

equation  (6.42)  on  page  182  from  Shapiro,  a relationship  between  the 
pressures  on  the  upstream  and  downstream  side  of  the  vent  can  be  developed. 
Doing  this  we  get: 


P,2-P2  = It^I  RT  1 2 (p./p) 


+ 4f§ 


(13) 


By  solving  equation  (13)  for  Mv,  differentiating  equation  (12)  and 
substituting  equations  (12)  and  (13)  into  equation  (9) , and  then  dividing 
both  sides  of  the  resulting  equation  by  Q/V  the  following  expression  is 


obtained: 


V.  + Qt 
l 


d (p/Pa)  _ _ £ 


d(t/t  ) 
e 


[i  - (p/pa)2l 

(Q/ A)  j4fL/D+2  ln(p  /pp 


(14) 


Equation  (14)  is  subject  to  the  initial  condition  p/p  = 1 at  t = 0 

ft 

and  can  be  numerically  integrated.  The  key  parameters  of  equation  (14) 

are  Q/A,  4f  L/D,  t/t  and  V. . A useful  and  more  easily  used  solution  to 

e i 

equation  (14)  can  be  arrived  at  by  assuming  d(p/p  )/d(t/t  ) * 0.  The 

cl  0 

rationale  for  equating  the  left  side  of  equation  (14)  to  zero  is  based  on 
ejqserimental  results,  where  it  was  observed  that  the  vacuum  magnitude 
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in  the  tank  would  approach  a finite  value  (depending  on  punping  rate  and 

vent  restriction)  and  remain  constant  within  the  accuracy  of  the  Data 

Acquisition  System  (DAS),  e.g,,  d(p/p  )/d(t/t  )=0.  Thus  eqt  ..ting  the 

a e ' 

left  side  of  equation  (14)  to  zero  results  in: 


£_  = 


RT 


RT  + (Q/ A)2  [2  In (p  /p)  + 4fL/D) 


(15) 


Equation  (15)  cannot  be  solved  directly,  but  must  be  iterated  to 


arrive  at  a solution.  Solving  equation  (15)  for  p/p  without  the  term 

a 

2 ln(p  /p)  in  the  denominator  and  then  using  this  first  step  value  of 
& 

p/p  in  the  complete  expression  and  resolving  for  p/p  , hastens  the 

d 3. 

iterative  process.  Convergence  is  rapid  since  2 ln(p  /p)  is  small  com- 

SL. 

pared  to  4f  L/D.  Thus,  equation  (15)  presents  a simplified  expression 
for  arriving  at  the  maximum  vacuum  to  be  expected  for  any  combination 
of  loading  rate  and  vent  restriction.  Additionally,  this  expression 
provides  a conservative  worst  case  for  cargo  offloading  in  that  cargo 
evaporation,  which  would  tend  to  decrease  the  maximum  tank  vacuum  attained, 
was  neglected  in  the  problem  formulation. 

B.  Experimental  Results  Unloading  Cargo 

In  order  to  compare  the  pressure-time  histories  for  various  vent 
systems  and  unloading  rates,  the  following  cargo  transfer  parameters  were 
chosen: 


4f  L/D  = 10,  20,  and  35 

Q/A  = 30,  50,  100  and  150  ft/sec. 

As  mentioned  previously  in  the  analytical  formulation,  the  effect  of 
cargo  evaporation  was  neglected,  (K=0)< providing  a conservative  worst 
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case  since  any  cargo  evaporation  will  compensate  for  a portion  of  the  vacuum 
buildup.  Additionally,  the  unloading  evaluations  were  performed  with  the 
tank  being  partially  full  (V^=0,  1/2,  1/3,  etc.)  to  assess  the  effect  of 
this  parameter  on  the  rate  of  vacuum  rise. 

Figure  14  shows  the  effect  the  initial  amount  of  cargo  in  the  tank 
has  on  the  rate  of  vacuum  rise.  From  Figure  14  it  can  be  seen  that  the 
smaller  the  vapor  space  is  when  offloading  is  begun,  the  more  rapid  is  the 
vacuum  buildup,  and  experimental  results  also  indicate  that  the  initial 
cargo  volume  does  not  affect  the  maximum  vacuum  attained  (within  the 
accuracy  of  the  DAS  and  for  the  range  of  parameters  explored) . Table  4 
presents  the  results  obtained  from  the  experimental  evaluation  of  cargo 
unloading  in  the  form  of  the  maximum  vacuum  attained  for  the  Q/A  and  4fL/D 
values  examined.  Also  included  in  Table  4 are  the  maximum  vacuums  as  pre- 
dicted by  the  simplified  analytical  model,  equation  (15).  One  will  note 
that  the  correlation  between  the  analytical  and  experimental  results  is 
quite  good  and  most  certainly  within  the  bounds  of  experimental  error. 

Figure  15  shows  the  relationship  between  cargo  unloading  rate  (.Q/A) 
and  maximum  vacuum  attained  in  the  tank  for  various  vent  system  restric- 
tions (4f  L/D)  as  predicted  by  equation  (15).  It  was  noted  during  the 
ejqjeriments  that  a buckling  of  the  tank  top  can  occur  as  a result  of 
the  combined  effects  of  the  tank  vacuum  created  by  unloading  cargo  and 
the  .geometry  of  the  tank  top, especially  when  the  tank  top  is  cambered. 

This  instability  may  occur  at  a vacuum  which  is  less  (in  absolute 
magnitude)  than  the  positive  pressure  required  to  cause  a yielding  type 
tank  failure.  The  combinations  of  tank  internal  vacuum  and  tank  top 
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FIGURE  14  - EFFECT  OF  INITIAL  LIQUID  VOLUME  OK  PRESSURE  - 
TIME  HISTORY  UNLOADING  CARGO  FOR  EXPERIMENTAL 
MODEL 
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TABLE  4 

ANALYTICAL  AND  EXPERIMENTAL 
MAXIMUM  VACUUM  UNLOADING  CARGO 


Q 

A 

(ft/sec) 

30 

Experimental 

Analytical 

50 

Experimental 

Analytical 

100 

Experimental 

Analytical 

150 

Experimental 

Analytical 

‘ACTUAL 


(12.6) 


NOMINAL  4ft 


‘(ACTUAL 


(25.2) 


Experimental  -.23  ^11-4^  -.44  ^22‘8^  -.77  ^41'0^ 
Analytical  -.25  -.45  -.78 


-1.38 

-1.35 


(19.5) 


‘(ACTUAL) 

35 

-.31 

(44.2) 

-.30 

-.77 

-.78 

(41.0) 

-2.19 

-2.18 

(34.0) 

(9.2) 


Actual  value  used  ejQ>erimentally  and  for  analytical 
calculations  (see  Appendix  A). 

Unloading  rate  and  vent  restriction  combination  produce 
choked  flow. 


geometric  configuration  which  would  cause  the  tank  and/or  deck  to  collapse 
inward  are  beyond  the  oonfines  of  this  investigation,  but  should  be  con- 
sidered when  establishing  safe  unloading  rates,  especially  for  tanks  with 
camber. 


Equation  (4)  shows  the  rate  of  tank  top  pressure  rise  during  liquid 


overfill  to  be  a function  of  the  following  parameters: 

1 dV 

B = tank  volume  of  expansion  coefficient  by  pressure  (—  — ) 

t p 

V^/A  = ratio  of  tank  volume  to  vent  area 

Q/A  = ratio  of  loading  rate  to  vent  area  > 

V-  = compressibility  of  liquid 

U = mean  velocity  of  liquid  in  the  vent  pipe 

The  volume  of  expansion  coefficient  for  the  model  tank  is  0. 0003B8 

psi  1 which  is  within  the  range  of  values  for  typical  prototype  tanks 

-1  1 

0.001  to  0.00001  psi  ) . The  compressibility  of  the  liquid  being  pumped 
into  the  tank  (water)  is  0.0000033  psi  1 and  therefore  is  generally  less 
than  the  tank  bulk  modulus;  hence  the  accomodation  of  pressure  rise  in 
the  tank  is  usually  governed  by  the  tank  expansion.  The  relationship 
between  the  model  cargo  compressibility  and  the  model  tank  bulk  modulus 
is  not  unlike  that  to  be  expected  from  a prototype  situation,  in  that  for 


the  case  of  a full  scale  tank  being  overfilled,  the  rate  of  pressure  rise 


would  be  governed  by  the  tank  bulk  modulus  more  so  than  the  compressibility 
of  the  cargo  (B»3f) . The  loading  rate  was  held  oonstant  during  each  test, 
and  loading  rates  (Q/A)  of  30,  50,  and  100  ft/sec  were  examined.  The 
ratio  of  tank  volume  to  vent  area  (V^A)  for  the  experimental  model  is 
40,196  ft  and  is  in  the  range  of  typical  V^/A  tank  volume  to  vent  area  ratios 
for  most  barges  (V^/ A barges  = 14,000  to  224,000  ft).  Thus,  the  overfill 
model  represents  a barge  tank  during  liquid  overfill  and  has  a pressure- 


time  history  which  corresponds  to  a barge  with  the  same  parameters.  Fig- 
ure 16'  shows  the  experimental  results  of  the  overfill  cases  run  for  load- 
ing rates  of  30,  50,  and  100  ft/ sec.  Also  included  in  the  figure  is  the 
expected  regime  of  tank  failure  for  vessels  of  this  type.  For  tins  case, 
tank  failure  can  be  expected  within  approximately  10  seconds  from  the 
t^me  the  tank  becomes  liquid  full.  The  tank  flexibility  plays  an  impor- 
tant role  in  helping  to  delay  the  pressure  buildup.  This  phenomenon  is 
due  to  the  fact  that  initially  the  outflow  velocity  of  the  liquid  is 
appreciably  less  than  the  inflow  velocity,  and  considerable  liquid 
accumulates  in  the  tank  due  to  tank  flexing  (e.g.  an  infinitely  stiff 
tank,  B=*0,  would  result  in  an  infinitely  rapid  pressure  rise),  since 
the  pressure-time  histories  generated  analytically  and  experimentally 
had  different  parameter  values,  a direct  comparison  of  results  could  not 
be  made.  However,  both  experimental  and  analytical  models  indicate  that 
tank  failure  will  occur  a short  time  after  the  tank  becomes  liquid  full 
(generally  less  than  one  minute). 

Equation  (6)  is  a steady  incompressible  flow  formulation  of  flow 
through  the  vent  system  and  defines  the  liquid  efflux  velocity  at  which 
cargo  can  pass  indefinitely  through  the  vent  system.  For  a typical  vent 
system  and  tank  configuration  with  4f  L/D=10,  a vent  pipe  rise,  H = 8ft 
and  an  allowable  pressure  P =6  psig,  the  maximum  loading  rate  to  vent 

ITlclX 

cross  sectional  area  is  about  6 ft/sec.  Thus,  for  loading  rates  higher 
than  this  (typically  higher  rates  are  employed),  tank  failure  will  occur 
in  a relatively  short  time  if  overfill  is  allowed  to  occur.  Because  of 
this,  some  type  of  device  to  protect  the  tank  should  be  employed. 


45 


B . Pre c -fc  His  tor1/  When  Vent  is  Blocked 

The  effect  of  a blocked  vent  or  stuck  P-V  valve  on  the  tank  pressure 
buildup  can  best  bo  examined  by  looking  at  two  separate  cases.  The  first 
would  be  when  the  vent  becomes  blocked  after  the  tan.,  is  liquid  full  and 
the  second  would  be  to  examine  the  effect  of  vent  blockage  before  the  tank 
is  full.  The  blocked  vent  was  modeled  by  installing  a ball  valve  in  the 
model  -rent  and  closing  it  off  at  the  appropriate  time  (e.g.  either  before 
or  after  overfill). 

For  the  case  of  a blocked  vent  after  overfill  had  commenced,  the  tank 
was  allowed  to  overfill  and  reach  a steady  state  flow  condition  (a  large 
diameter  vent  pipe  was  employe  * and  the  resulting  steady  state  pressure 
was  less  than  1 psig)  and  the  ball  valve  was  than  closed,  with  time 
beginning  after  the  valve  was  closed.  Even  though  the  tank  was  allowed 
to  become  liquid  full,  a small  amount  of  trapped  air  remained  in  the 
tank  due  to  the  tank  top  geometry  and  the  intrusion  of  the  vent  pipe 
fixture  back  into  the  tank.  The  effect  of  the  trapped  air  on  the  rate  of 
pressure  rise  in  the  tank  was  examined  by  running  the  blocked  vent  cases 
with  the  tank  physically  oriented  in  two  ways  as  is  depicted  in  Figure  17. 

The  effect  that  the  amount  of  trapped  air  had  on  the  rate  of  pressure 
rise  in  the  tank  can  be  observed  in  Figure  18  where  the  experimental 
pressure  time  histories  for  both  barrel  orientations  are  shown  along 
with  the  analytical  pressure  history  for  the  same  case  with  no  trapped 
air.  This  discrepancy  between  the  results  is  due  to  the  entrapped  air  in 
the  tank  and  can  be  best  explained  by  considering  tne  stiffness  of  the 
tank  and  the  compressibility  of  the  trapped  air  as  two  springs  in  series. 
In  this  case  though,  one  spring  stiffness  (the  air)  is  initially  mucn 
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TILTED  CONFIGURATION 
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UPRIGHT  CONFIGURATION 

FIGURE  17  - EFFECT  OF  TANK  ORIENTATION  ON  THE 
AMOUNT  OF  TRAPPED  AIR  IN  THE  TANK 


(PSIG) 


softer  and  has  a stiffness  which  is  changing  rapidly  (becoming  stiffer) 
with  respect  to  the  other  (the  tank)  as  the  pressure  inside  the 
tank  is  increasing.  The  relationship  between  stiffness  and  bulk  modulus 
is  that  a stiffer  object  would  have  a smaller  bulk  modulus.  When  the 
stiffness  of  the  air  becomes  greater  than  that  of  the  tank  structure 
itself,  the  stiffness  of  the  tank  would  then  begin  to  govern  the  pressure 
rise  in  the  tank  and  approach  the  rate  of  pressure  rise  predicted  by  equa- 
tion (5)  which  does  not  consider  air  trapped  in  the  tank. 

Than,  by  running  the  blocked  vent  case  with  the  tank  in  a tilted 
configuration  (little  trapped  air) , a pressure-time  history  is  attained 
which  more  closely  resembles  that  predicted  by  the  analytical  model  than 
that  obtained  by  running  the  test  with  the  tank  in  a upright  configu  \- 
tion.  That  is,  the  slope  of  the  tilted  configuration  p vs.  t curve  more 
rapidly  approaches  the  slope  of  the  analytical  curve  than  does  the  p vs.  t 
curve  of  the  vertical  configuration.  Further,  it  is  felt  that  both  curves 
(tilted  and  vertical  configurations)  would  have  eventually  attained  the 
slope  of  the  anlytical  model  had  the  test  not  been  terminated  at  25  psig 
for  reasons  of  safety. 

Figures  19  and  20  show  experimental  pressure-time  histories  for  two 
vent  restrictions  (4f  L/D  = 10, oo  )and  two  loading  rates  (Q/A  = 100  and  23 
ft/sec)  and  analytical  pressure-time  histories  for  the  infinite  vent  restric- 
tion or  blocked  vent  case.  One  will  note  from  Figures  19  and  20  that  the 
difference  in  the  pressure-time  history  fot  a finite  (4f  L/D  - 10)  and 
infinite  vent  restriction  (blocked  vent)  is  relatively  small  in  the  early 
stages  of  overfill  up  to  the  regime  of  tank  failure.  Thus,  the  blocked 


50 


TANK  TOP  PRESSURE  (PSIG) 


vent  case  provides  an  easily  employed,  although  conservative,  worst  case 
solution  for  estimating  the  time  to  tank  failure  once  overfill  begins. 

The  second  case  to  consider  is  that  of  a vent  being  blocked  before 
the  tank  becomes  liquid  full  (e.g.  a stuck  PV  valve).  This  situation  was 
modeled  experimentally  by  closing  the  ball  valve  in  the  vent  system  before 
the  tank  was  liquid  full.  Figure  21  shows  an  experimentally  derived 
pressure-time  history  arrived  at  in  the  manner  just  described.  The  ball 
valve  was  closed  with  the  tank  60%  full  (simulating  a relief  valve  incap- 
able of  relieving  tank  pressure)  employing  a simulated  transfer  rate 
Q/A  of  23  ft/sec. 

The  phenomena  occurring  is  the  same  as  that  described  previously  when 
there  was  air  trapped  in  the  tank  during  liquid  overfill,  only  now  because 
of  the  large  volume  of  air  trapped  in  the  tank,  the  bulk  modulus  of  the 
trapped  air  is  net  changing  as  rapidly  as  before  and  thus  the  pressure  rise 
is  not  as  rapid.  As  one  can  see  from  Figure  21,  where  the  tank  pressure 
is  plotted  aginst  the  fraction  of  total  fill  time,  the  pressure  rise  in  the 
tank  approaches  the  failure  regime  of  the  tank  long  before  the  tank  becomes 
liquid  full.  Thus,  the  employment  of  warning  devices  based  on  the  tank 
liquid  level  would  be  ineffective  against  this  situation.  An  automatic 
shutoff  mechanism  triggered  by  pressure  would  be  needed  to  provide  protec- 


tion against  this  occurrence. 


TANK  TOP  PRESSURE  (PSIG) 


25 


FIGURE  21 


EFFECT  OF  BLOCKED  VENT  BEFORE  TANK  BECOMES  LIQUID  FULL 
FOR  EXPERIMENTAL  MODEL  ^ LL 
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VII.  SUMMARY  OF  FINDINGS 


A.  Normal  Transfer  of  Cargo 

1.  The  analytical  model  of  cargo  loading  with  and  without  cargo 
evaporation,  equation  (2)  , as  developed  by  Wilson  and  Raj  ^"predicts 
the  pressure-time  history  of  a cargo  tank  being  loaded  within  the 
constraints  of  the  assumptions  made  by  the  analytical  model  for  evapora- 
tion . 

2.  An  analytical  model  for  cargo  offloading,  equation  (14),  can  be 
employed  to  predict  pressure-time  histories  for  a cargo  tank  being  off- 
loaded. This  formulation  is  conservative  for  volatile  cargoes.  A simpli- 
fied expression,  equation  (15) , was  derived  from  equation  (14)  to  predict 
the  maximum  vacuum,  attained  offloading  cargo,  and  has  been  validated 
through  scale  model  tests. 

3.  From  the  analyses  and  experimental  results,  the  following  conclu- 
sions were  reached  for  non-volatile  cargoes: 

• The  maximum  tank  pressure  for  loading  non-volatile  'cargo  can  be 
calculated  using  the  simplified  expression  for  loading  cargo,  equation  (8). 

• The  maximum  tank  vacuum  occurring  when  offloading  cargo  can  be 
calculated  using  the  simplified  expression  for  offloading  cargo,  equation 
(15). 

• The  maximum  tank  pressures  resulting  during  the  loading  of  non- 
volatile cargoes  are  small  compared  to  calculated  lank  failure  pressures 
for  typical  loading  rates  and  vent  restrictions  (Q/A  = 30  ft/sec  and  4f 
I/D  = 10)  . 

• the  absolute  value  of  the  maximum  vacuum  attained  unloading  cargo 
is  less  than  the  absolute  value  of  the  maximum  tank  top  pressure  attained 
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when  loading  cargo,  indicating  the  danger  of  tank  failure  defined  as 


reaching  tank  material  yield  to  be  less  than  when  loading  cargo,  but 
the  possibility  of  tank  top  buckling  at  a lower  pressure  (due  to  tank 
top  geometry  conforming  to  that  of  a cambered  deck)  when  unloading 
cargo  should  also  be  considered. 

• For  high  loading  rates  and  vent  restrictions  the  maximum  tank  top 
pressure  changes  rapidly  with  small  changes  in  loading  rate,  indicating 
the  danger  in  employing  these  higher  rates. 

4.  For  normal  vent  systems,  cargo  volatilities,  and  cargo  loading 
rates,  it  appears  that  volatile  cargoes  can  be  loaded  with  an  adequate 
safety  factor  for  tajik  pressure  rise  (subject  to  the  assump.  -ons  of  the 
evaporation  model) . 

B.  Liquid  Overfill 

The  steady  incompressible  flow  formulation  for  flow  through  the  vent 
system  as  given  in  equation  (6)  r ^ 


Q/A  = U = 


p-pa7*gH 


h yopf  L/EH-gj 

can  be  used  to  determine  the  maximum  allowable  loading  rate  to  vent  cross 
section  area  ratio  (2/A  = U)  at  which  cargo  could  safely  pass  indefinitely 
through  the  vent  system.  For  a tank  failure  pressure  of  6 psig  with  a vent 
restriction  of  10  and  a vertical  rise  of  8 ft,  the  maximum  Q/A  ratio 
would  be  approximately  6 ft/sec.  For  a higher  rate  of  loading  or  a smaller 
vent  system  capacity,  tank  failure  will  occur  within  a relatively  short 
time  after  the  tank  becomes  liquid  full. 

The  effect  of  trapped  air  in  the  tank  has  the  net  effect  of  decreasing 
the  rate  of  pressure  rise  in  the  tank,  but  not  to  the  extent  of  diminishing 
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the  danger  presented  by  liquid  overfill.  Tank  flexure  plays  a similar 
role  in  relieving  overfill  pressure,  with  a stiffer  tank  reaching 
failure  pressures  more  rapidly  than  a more  flexible  one.  Further,  with 
these  considerations,  it  still  appears  that  tank  failure  will  occur 
within  about  30  seconds  after  the  tank  becomes  liquid  full,  for  normal 
combinations  of  loading  rate  and  vent  restriction. 

Thus,  it  appears  that  some  type  of  automatic  shut-off  mechanism 
triggered  by  pressure  to  shut  off  the  liquid  flowing  into  the  tank  must 
be  provided  to  protect  the  tank.  Another  possibility  would  be  a 
pressure  triggered  hatch  cut  into  the  tank  top  at  the  deck  level  which 
has  sufficient  area  to  allow  liquid  overflow  without  any  dangerous 
pressure  buildup.  In  addition,  any  system  employed  should  have  high 
level  alarms  to  warn  the  operator  of  the  impending  danger. 


VIII.  APPLICATION  TO  VENT  SYSTEM  EVALUATION  AND  DESIGN 


One  can  evaluate  the  adequacy  of  an  existing  vent  system  for  avoid- 
ing overpressure  in  an  individual  tank  by  errploying  the  following  steps: 

1.  Calculate  the  characteristic  venting  velocity  Q/A  (ft/sec)  for 

• 3 

the  tank  in  question  by  determining  the  loading  rate  Q (ft  /sec)  from 

the  loading  rate  and  number  of  tanks  to  be  loaded  simultaneously  and 

2 

dividing  by  the  cross-sectional  area  A (ft  ) of  the  vent  pipe (s)  coming 
out  of  the  tanks. 

2.  Calculate  the  effective  length  to  diameter  ratio  1/D  of  the  vent 
system  by  first  dividing  each  pipe  length  employed  in  the  vent  system  by 
its  respective  diameter  and  then  summing.  Then  add  in  all  the  equivalent 
L/D  ratios  for  all  bends,  elbows,  tees,  valves,  flame  arrestors  and  the 
entrance  loss  associated  with  each  vent  pipe.  These  values  can  be  found 
in  most  handbooks  or  References  1,  5 or  6. 

3.  The  maximum  expected  tank  pressure  loading  cargo  can  be  determined 
using  Table  5 and  the  values  of  Q/A  and  L/D  determined  in  steps  1 and  2. 
Table  5 was  developed  by  ADL  using  equation  (2)  with  a high  evaporation 
rate  (K  = 0.7),  which  it  is  felt  provides  some  factor  of  safety. 

4.  The  maximum  expected  vacuum  created  in  the  tank  when  offloading 

cargo  can  be  determined  using  equation  (15)  and  the  values  of  Q/A  and  L/D 

2 2 

determined  in  steps  (1)  and  (2)  (assume  4f  = .02  and  RT  = 900,000  ft  /sec  ). 

5.  The  maximum  pressure  during  cargo  overfill  can  be  arrived  at  by 
.first  determining  the  height  of  the  vertical  portion  of  the  vent  pipe  H 
and  then  using  equation  (6)  and  the  values  of  H,  Q/A  and  4f  L/D  previously 
determined. 
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-ABLE  5 

MAXIMUM  TANK  PRESSURES  (PSIG)  for 
GIVEN  VALlTto  OF  VENT  SYSTEM 
L/D  and  Q/A  (BASED  ON  f = .005) 
K •’  0.7 

(FROM  R ^ERENCE  1) 


Q/A  (FT/SEC) 


600 

400 

1000 

1200 

1600 

2000 

24(10 

3200 

4()U0 

6000 

WJ'J 

10000 

12000 

16000 

2 ijV) 

24000 

3 2000 


O'- 


Tank  pressure*  less  than 
3 psi  y 


2.2  3.0 

2.5  3.5 

2.7  3.8 

2.9  A. 2 

3.3  4.8 

3.6  5.4 

3.9  6.0 
4.5 


3.3  4.0  4.8 
4.0  4.9  6.0 
4.6  5.8 


3.0  3 . 


/ 2.8  3.6 

/ 3.1  4.1 

3.0  3.5  4.6 

3.2  3.8  5.0 

3.6  4.3  5.9 


2.7 

3.2 

4.2 

3.2 

3.8 

5.1 

3. 

4.4 

5.9 

4.1 

4.9 

4.4 

5.4 

/ 

5.1 

/ 

5.8 

/ 

Tank  pressure  greater 
than  6 psi  - hazardous 
for  most  vessels 


6.  Compare  the  expected  pressures  arrived  at  in  steps  3 through  5 
with  the  expected  tank  failure  pressure.  Tank  failure  pressure  can  be 
obtained  from  a structural  analysis  of  the  tank  in  question  or  in  lieo 
of  such  an  anlysis  the  ADL  value  of  6 psig  may  be  employed. 

Based  on  estimates  of  tank  failure  pressure  and  typical  vent  re- 
strictions ( 6psig  and  4f  h/D  = 3 to  10) , experimental  results  indicate  that 
typical  vent  systems  in  use  today  have  an  adequate  capacity  for  gas  vent- 
ing but  an  inadequate  capacity  for  liquid  overfill.  Further,  any  vent 
system  designed  to  facilitate  liquid  overfill  v^ould  be  extremely  overde- 
signed for  normal  gas  verting.  Thus,  it  appears  that  a practical  approach 
to  vent  system  design  would  be  to  design  the  vent  system  to  accomodate 
normal  gas  venting  with  an  adequate  safety  factor  in  addition  to  incorpora- 
ting preventative  measures  in  the  tank/vent  system  to  avoid  or  relieve  the 
occurence  of  liquid  overfill. 
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NOMENCLATURE 


Units 


A 


B 


f 

9 

H 

K 


L 


L/D 


M 

v 

M 

c 

M 

vap 


M 

M. 

X 

p 

Pa 

P 

6 

R 


Re 

t 


fill 


=■  cross  sectional  flow  area 

= coefficient  of  volume  expansion  of  tank  by  pressure 


= 1-  (dV) 

VT 

= fanning  friction  factor  in  the  vent  pipe 
= acceleration  due  to  gravity 

= height  of  the  vertical  portion  of  the  vent  pipe 

= ratio  of  evaporation  rate  to  the  loading  rate 
taken  at  tfj_n 


= effective  length  of  vent  pipe 

= effective  length  to  diameter  ratio 

= rate  of  mass  efflux  through  the  vent 

= characteristic  mass  efflux  rate  of  air  =y3^  Q 

= rate  of  vapor  mass  addition  to  the  ullage  by 
evaporation 

= mass  of  liquid  in  the  tank  at  any  time 
= rate  of  mass  inflow  (liquid  filling  rate) 

= pressure  inside  the  tank  (at  the  top  of  the  tank) 
= atmospheric  pressure 
= gauge  pressure  inside  the  tank 
= volumetric  liquid  filling  rate 
= gas  constant  - 1698.7 

= Reynolds  number  for  flow  in  the  pipe  = ~ 

= time 

= characteristic  tank  filling  time  (volume  of  tank 
divided  by  loading  rate) 


ft 


.-1 

psi 


2 

ft/sec 

ft 

ft 

lbm/sec 

lbm/sec 

lbm/sec 

lbm 

lbm/sec 

psi 

psi 

psig 

f t 2/sec 

ft2/sec2/OR 

sec 

sec 
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NOMENCLATURE  (CONT) 


Units 


T 

U 

V 


V. 

1 

VT 

c 

X 

/T 

a 

A 

/vap 


= temperature 

= mean  velocity  of  liquid  in  the  vent  pipe 

= tank  volume  not  occupied  by  liquid  at  any  instant 
of  time 

= initial  tank  volume  not  occupied  by  liquid 
= total  tank  volume 
= vent  pipe  roughness 
= compressibility  of  liquid  = ^ ^ 

= kinematic  viscosity 
= density  of  liquid 
= density  of  gas 
= density  of  pure  vapor 


ft/sec 

ft3 


f li- 
ft3 
ft 

.-1 

psi 

2 

ft  /sec 
lbm/ft3 
lbm/ft3 


lbm/ft 


6 


REFERENCES 


Wilson,  R.  P. , Jr.  and  Raj/P.  K.  P.,  Vent  System  and  Loading  Criteria 
For  Avoiding  Tank  Overpressurization,  Department  of  Transportation 
Report  No.  CG-D-59-77,  September  1977. 

Hovanesian,  J.  and  Y.  Y.  Hung,  Moire'’  Contour-Sum  Contour-Difference 
and  Vibration  Analysis  of  Arbitrary  Objects,  Applied  Optics , Vol.  10, 
No.  12,  December  1971. 

Meadows,  D.  M. , W.  0.  Johnson  and  J.  B.  Allen,  Generation  of 
Surface  Contours  by  Moire  Patterns,  Applied  Optics , Vo.  9,  No.  4, 
April  1970. 

Shapiro,  A.  H. , The  Dynamics  and  Thermodynamics  of  Compressible 
Fluid  Flow,  Vol.  I,  P.  83,  Ronald  Press,  New  York  (1953). 

Baumeister,  (Editor) , Marks'  Handbook  for  Mechanical  Engineers, 

7th  Edition,  McGraw-Hill,  New  York,  (1967). 

Flow  of  7 1uids  Through  Valves,  Fittings,  and  Pipe,  Crane  Co. , 
Technical  Paper  No.  410,  Chicago  (1969). 

Streeter,  V.  L. , Fluid  Mechanics , McGraw-Hill,  4th  Edition, New 
York  (1966). 


1 


i 


* - 

i • 

APPENDIX  A 

DETERMINATION  OF  MODEL  VENT  PIPE  ROUGHNESS 
Accurate'  determination  of  the  model  vent  pipe  roughness  was  necessary 
to  accurately  model  the  vent  system  restriction,  4f  L/D.  Since  a model 
vent  pipe  of  diameter  D could  be  cut  to  any  length  L,  accurate  modeling  of 
the  vent  restriction  would  be  directly  related  to  the  establishment  of  the 
pipe  Fanning  friction  factor  f.  The  friction  factor  for  a particular 
pipe  is  a direct  function  of  Reynolds  number  Re,  and  pipe  relative  roughness 
<-/D,  where  the  two  parameters  can  be  related  on  a Moody  diagram.  In  a 
typical  pipe  flow  calculation  all  quantities  except,  4f,  can  be  measured. 

Thus,  by  conducting  a pipe  flow  experiment  using  the  model  vent  pipe,  one 
would  be  able  to  calculate  the  friction  factor ,4f , directly  for  the  par- 

} 

ticular  Reynolds  number  at  which  the  test  was  run.  Knowing  these  two 
quantities  and  employing  a Moody  diagram,  one  would  then  be  able  to  de- 
termine the  model  vent  pipe  relative  roughness  and  enable  one  to  model  the 
vent  restriction  4f  L/D  for  any  Reynolds  number  (e.g.  gas  efflux  velocity 
through  the  vent  system) . 

An  attempt  was  made  to  arrive  at  the  vent  pipe  roughness  by  measuring 
the  head  loss  resulting  with  water  flowing  through  the  model  vent  pipe. 

But  because  of  the  inability  to  attain  a steady  state  of  flow  and  the  lack 
of  consistency  of  results,  this  method  of  roughness  determination  was 
abandoned  in  favor  of  the  following  method.  A known  length  of  model  vent 
pipe  was  installed  in  the  tank  top  and  air  was  then  pumped  through  a flow- 
meter at  a constant  rate  venting  to  the  ambient.  The  pressure  on  the 
upstream  side  of  the  vent  pipe  was  monitored  by  a pressure  transducer 
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located  on  the  tank  top.  A sketch  of  the  pipe  roughness  experimental 
test  setup  can  be  found  in  Figure  A-l. 

Sy  writing  an  energy  balance  between  points  1 and  2 and  accounting 
for  all  losses  the  following  equation  results: 


4f  =» 


KHI'S-’M’-'.) 


D_ 

L 


(Al) 


where  C is  the  orifice  contraction  coefficient,  accounting  for  the  losses 
c 

of  the  connecting  fixture  at  the  vent  pipe  orifice.  Equation  (A  1)  con- 
tains two  unknowns,  4f  and  C . By  running  the  experiment  with  tvro  lengths 

c 

of  pipe,  the  resulting  equations  could  be  solved  simultaneously  for  4f  and 


C . 
c 

Results  of  the  roughness  evaluations  run  on  the  pipe  sections  to  be 
used  in  the  venting  model  are  shown  in  Figure  A-2  as  points  on  a Moody 
diagram.  Taking  the  average  of  the  relative  roughness  values  from  these 
evaluations,  yielded  a relative  roughness  C/D  =0.0025  or  a pipe  roughness 
£ = .000038  ft.  This  roughness  value  lies  between  that  of  commercial 
steel  and  drawn  tubing  (C  = .00015  and  €=  .000005  respectively).  The 
model  vent  pipe  material  is  drawn  stainless  steel. 

This  relative  roughness  value  was  used  to  model  the  vent  restrictions 
employed  for  subsequent  evaluations.  The  vent  restriction  was  modeled  by 
first  establishing  the  Reynolds  number  at  which  the  test  would  be  run. 

This  would  simply  be  tile  vent  pipe  velocity  of  the  test  (Q/A)  times  the 
model  vent  pipe  diameter  (D)  divided  by  the  kinematic  viscosity  of  the 
fluid  being  vented  (e.g.  air).  The  friction  factor,  4f,would  then  be  es- 
tablished by  finding  the  intersection  of  the  relative  roughness  curve  for 
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FIGURE  A1  - SCHEMATIC  OF  EXPERIMENTAL  TEST  SETUP 

FOR  VQJT  PIPE  ROUGHNESS  DETERMINATIONS 
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the  pipe  and  the  Reynolds  number  at  which  the  test  was  run  on  a Moody 
diagram.  Knowing  the  friction  factor  4 f and  the  vent  pipe  diameter  D, 
the  desired  value  of  4f  I/D  could  now  be  modeled  (for  this  particular 
Reynolds  number)  by  cutting  tne  model  vent  pipe  to  the  appropriate 
length,  L. 

For  example: 

To  model 

Q/A  = 100  ft/sec  and  4f  L/D  = 20  with  D = .0155  ft  compute. 

Re  = VD /ff  = 100  x .0155/.  000162  = 9568,  from  Moody  diagram  for 
Re  * 9568  and  </D  » .0025  we  get 
4f  = .035  and  then  L = 8.85  ft. 

Since  the  flow  state  for  the  evaluations  was  not  fully  turbulent 
(4f  not  constant)  and  new  lengths  of  pipe  were  not  employed  for  each 
loading  rate  examined,  the  actual  vent  restriction  4f  L/D  was  arrived  at 
using  the  pipe  relative  roughness  to  account  for  4f  not  being  constant 
with  changing  loading  rates.  Consider  the  above  example,  but  with  the 
vent  pipe  velocity  Q/A  = 50  ft/sec.  For  this  case  the  Reynolds  number 
Re  = 4784  and  the  friction  factor  4f  = .041,  resulting  in  4f  L/D  = 22.8. 

This  correction  is  reflected  in  the  (actual)  4f  L/D  values  given  in 
Tables  2,  3 and  4.  Additionally,  in  the  design  of  the  transfer  experimental 
parameters,  care  had  to  be  taken  so  that  a combination  of  4f  L/D  (pipe 
length  ar.d  diameter)  and  transfer  rate  Q/A  (punping  rate  and  vent  area) 
were  not  employed  such  that  Mach  1 was  attained  in  the  vent. 
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APPENDIX  B 


«OIRE/  METHOD  FOR  MEASURING  PLATE  DEFLECTIONS 

Since  95%  of  the  tank  volume  change  by  pressure  occurred  from  the 
tank  top  and  bottom  deflections  and  since  the  change  in  tank  volume  with 
pressure  (bulk  modulus)  is  an  important  parameter  for  overfill,  it  was 
felt  that  a backup  method  (in  addition  to  the  linear  potentiometer)  for 
measuring  the  tank  top  and  bottom  deflections  would  be  advantageous.  To 
accomplish  this,  the  Moire'  Contour  Sum-Contour  Difference  Method  was  em- 
ployed to  measure  the  plate  deflections  in  the  tank  top  and  bottom. 

Briefly,  the  Moire'  method  enployed  here  consisted  of  the  observation 
of  contour-like  patterns  created  by  the  mode  interference  of  a grid  with 
its  shadow  cast  onto  a surface  (the  tank  top) . A photograph  of  the  pattern 
is  made  and  stored  and  the  tank  top  is  then  allowed  to  deform  by  increasing 
the  internal  tank  pressure.  A second  exposure  of  the  tank  top  is  then  made 
over  the  first  exposure  by  double  exposing  the  film.  The  result  of  this 
process  is  the  development  of  Moira"  interference  fringes  which  depict 
contour  differences  on  the  surface  photographed.  The  contour  elevations 
can  be  calculated  knowing  the  distances  of  the  projected  grid  and  camera 
from  the  surface,  the  angles  they  form  with  the  surface  in  question,  and 
the  grid  spacings. 

Figure  B-l  shows  the  deflections  measured  on  the  tank  top  and  bottom 
employing  the  Moire'  method  to  arrive  at  the  change  in  tank  volume  for 
a 25  psig  internal  pressure.  These  measurements  were  made  statically  by 
projecting  the  Moire'  grid  on  the  tank,  taking  an  exposure  of  the  pattern, 
pressurizing  the  tank  to  25  psig  and  taking  a second  exposure.  The  resulting 
photograph  recorded  the  Moire  fringes  thus  produced,  which  when  converted 
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FIGURE  B1  - RESULTS  OF  STATIC  TANK  BULK  MODULUS  DETERMINATION 
USING  MOIRE"  METHOD 
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to  deflections  produced  the  deflection  pattern  shown  in  Figure  B-l.  A 
micrometer  was  used  to  measure  the  deflection  of  the  center  of  the  tank 
top  and  bottom  plate  to  give  an  additional  data  point  for  this  static  tank 
calibration. 

To  arrive  at  the  tank  bulk  modulus  B = AV/Ap  VT,  additional  infor- 
mation from  the  tank  strain  channels  (hoop  strain  and  the  longitudinal 
cylindrical  extension  strains)  was  employed  in  conjunction  with  the  Moire^ 
results,  to  arrive  at  the  total  tank  volume  change.  The  results  of  this 
static  bulk  modulus  determination  can  be  found  in  Table  B-l. 

The  value  for  B thus  determined  was  then  used  in  the  analytical 
expression  for  a blocked  vent  to  predict  the  blocked  vent  pressure  time 
history.  Additionally,  the  Moire^  method  was  enployed  during  the  actual 
experimental  overfill  evaluations  to  get  the  time  dependent  volume  changes 
during  overfill.  This  was  accomplished  by  exposing  motion  picture  film 
of  the  grid  on  the  tank  top  before  testing,  rewinding  the  exposed  film, 
then  running  the  test  and  re-exposing  the  same  film  during  the  test.  The 
result  was  a time  dependent  contour  map  which  could  be  equated  to  a time 
dependent  volume  change  and  bulk  modulus.  The  change  in  bulk  modulus  as 
a function  of  time  (pressure)  was  less  than  20%  and  the  assumption  of  a 
constant  value  for  the  bulk  modulus  considered  valid  for  comparative 
purposes  with  the  analytical  model  (assumes  constant  B) . 

TABLE  B-l 


Tank  Pressure  = 

25  psig 

V 

Tank 

= 13,281  in 

Av„  , „ 

69.0  in3 

A VT 

- 1 in 

Tank  Top 

Long 

A v = 

57.0  inJ 

A V 

= 2 in3 

Tank  Bottom 

Hood 

AV 

129 

B = 

vtaP 

13,281x25 

. 000388  psi 
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